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Abstract. In this lecture the basic concepts used in accretion disk theory are reviewed, with 
emphasis on aspects relevant for X-Ray Binaries and Cataclysmic Variables. 
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1. Introduction 


Accretion disks are inferred to exist in objects of very different scales, ranging from 
the neutron-star-to-solar scale in Low Mass X-Ray Binaries (LMXB) and Cata¬ 
clysmic Variables (CV), the solar-to-AU scale in protostellar disks, and the AU to 
parsec scale for the disks in Active Galactic Nuclei (AGN). On the kiloparsec scale 
galaxies, especially spirals share some properties with accretion disks and can be 
regarded as the only ‘disks’ whose structure is, at present, directly observable (the 
differences are sufhciently large however, that I will exclude them from further dis¬ 
cussion in this lecture). 

Lacking direct (i.e. spatially resolved) observations of disks, theory has tried to 
provide models, with varying degrees of success. Since some basic questions about 
disks are still unsolved (in particular the source of viscosity), the theory of disks 
is still a fairly undivided held. Progress made by observations or modeling of a 
particular class of objects is likely to have direct impact for the understanding of 
other objects. 

An interesting connection exists of accretion disks with jets (such as the spectac¬ 
ular jets from AGN and protostars), and outflows (the ‘GO-outflows’ from protostars 
and possibly the ‘broad-line-regions’ in AGN). A magnetically driven wind from the 
inner regions of an accretion disk is currently to most promising model for these 
phenomema. See Blandford (1989), Kdnigl and Ruden (1993), Wardle and Kdnigl 
(1993), Emmering et al. (1992), Spruit (1993). 

In this lecture I concentrate on the more basic aspects of accretion disks. Impor¬ 
tant subjects like theories for the viscosity in disks are dicussed only briefly (section 
8), because of complexity and lack of conclusiveness of the current state of the art. 
Emphasis is on those aspects of accretion disk theory that connect to the observa¬ 
tions of LMXB and GV’s. 

Eor other reviews on the basics of accretion disks, see Pringle (1981), Verbunt 
(1982), Erank, King and Raine (1992), Treves et ah, (1988). 
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2. Accretion: general 

Gas falling into a point mass potential 

GM 

$ =- 

r 


from a distance rg to a distance r converts gravitational into kinetic energy, by an 
amount GM(l/r — 1/ro) or for simplicity, assuming that the starting distance is 
large, GM/r. If the gas is then brought to rest, for example at the surface of a star, 
the amount of energy e dissipated per unit mass is 

GM 

e = - (res 

r 

or, if it goes into a circular orbit at distance r: 


e 


1 GM 

2 r 


(orbit). 


The dissipated energy goes into internal energy of the gas and into radiation which 
escapes to inhnity (usually in the form of photons, but neutrino losses can also play 
a role). 


2.1. Adiabatic accretion 

Consider Rrst the case when radiation losses are neglected. This is adiabatic accre¬ 
tion. For an ideal gas with constant ratio of specific heats 7 , the internal energy per 
unit mass is 

( 7 ^' 

With the equation of state 

P = TZpT/fi (1) 

where TZ is the gas constant, p the mean weight per particle, we find the temperature 
of the gas after the dissipation has taken place (assuming that the gas goes into a 
circular orbit): 

( 2 ) 

where Tyir, the virial temperature is given by 


Tvi 


GMp 

TZr 


In an atmosphere with temperature near Tyir, the sound speed is close to the escape 
speed from the system, the hydrostatic pressure scale height is of the order of r, 
and such an atmosphere evaporates on a relatively short time scale in the form of a 
stellar wind. 

A simple example is spherical adiabatic accretion, first discussed by Bondi (1952). 
The main finding is that such accretion is possible only if 7 < 5/3. This makes sense: 
the larger 7 , the larger the temperature in the accreted gas (eq. 2 ), and beyond a 
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critical value the temperature is too high for the gas to stay bound in the potential. 
A classical situation where adiabatic and roughly spherical accretion takes place is 
the supernova implosion: when the central temperature becomes high enough for 
the radiation held to start desintegrating nuclei, 7 drops and the envelope collapses 
onto the forming neutron star via a nearly static accretion shock. Another case are 
Thorne-Zytkov objects (see Ergma, this volume), where 7 can drop to low values 
due to pair creation, initiating an adiabatic accretion onto the black hole. 

Adiabatic spherical accretion is fast, taking place on the dynamical or free fall 
time scale 

Td = rjv ^ = (3) 

When radiative loss becomes important, the accreting gas can stay cool irrespec¬ 
tive of the value of 7 , and Bondi’s critical value 7 = 5/3 plays no role. With losses, 
the temperatures of accretion disks are usually much lower than the virial temper¬ 
ature. The optical depth of the accreting flow increases with the accretion rate M. 
When the optical depth becomes large enough so that the photons are ‘trapped’ 
in the flow, the accretion just carries them in, together with the gas (Rees 1978, 
Begelman 1979). Adiabatic accretion therefore occurs above a certain critical rate 
Me. To compute this rate, consider photons with a mean free path I given by 

I = (kp)“\ 

where k is the opacity (cm^/g) and p the density. Over a length scale r these photons 
diffuse outward at a speed r’diir: 


I'd® ~ Ic/r. 

Assuming spherical symmetry the accretion velocity is, in terms of the accretion 
rate: 

M 

p 

The critical rate is obtained when the outward diffusion speed just balances the 
accretion speed, which yields 

Me = Airrc/ K. 

This happens to be equal to the Eddington rate discussed below. In most cases. Me 
marks the boundary between adiabatic and cooling accretion (for special circum¬ 
stances where the critical rate is not the Eddington rate, see Spruit et al. 1987). 
In a fully ionized gas with radiation pressure negligible, 7 = 5/3, so that Bondi’s 
condition for adiabatic accretion is only marginally satisRed. If radiation pressure 
dominates, 7 is of the order 4/3, and adiabatic spherical accretion is possible. 

2.2. The Eddington Limit 

Objects of high luminosity have a tendency to blow their atmospheres away due 
to the radiative force exerted when the outward traveling photons are scattered or 
absorbed. Consider a volume of gas on which a flux of photons is incident from one 
side. Per gram of matter, the gas presents a scattering (or absorbing) surface area of 
K cm^. The force exerted by the radiative flux F on one gram is Fk/c. The force of 
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gravity pulling back on this one gram of mass is GMjr’^. The critical flux at which 
the two forces balance is 


Fe 


c GM 


( 4 ) 


Assuming that the flux is spherically symmetric, this can be converted into a critical 
luminosity 

Le = 'iirGMc/K, (5) 


the Eddington luminosity (e.g. Rybicki and Lightman, 1979). If the gas is fully 
ionized, its opacity is dominated by electron scattering, and for solar composition 
K is then of the order 0.3 cm^/g (about a factor 2 lower for fully ionized helium, a 
factor up to 10^ higher for partially ionized gases). With these assumptions. 


Le 


M 


1.710-'''- erg/sf«410'‘ 

Mq 


M 


L(7) 


If this luminosity results from accretion, it corresponds to the Eddington accretion 
rate Me'- 

=Le Me = Mrc/n. (6) 

r 

Whereas Le is a true limit that can not be exceeded by a static radiating object 
except by geometrical factors of order unity (see chapter 10 in Erank et al, 1992), no 
maximum exists on the accretion rate. Eor M > Me the plasma is just swallowed 
whole, including the radiation energy in it (cf. discussion in the preceding section). 
With K = 0.3: 

ME~l-310'®r6 g/sRi210“®r6 Mgyr"', 
where re is the radius of the accreting object in units of 10 km. 


3. Accretion with Angular Momentum 

When the accreting gas has a zonzero angular momentum with respect to the accret¬ 
ing object, it can not accrete directly. A new time scale appears, the time scale for 
outward transport of angular momentum. Since this is in general much longer than 
the dynamical time scale, much of what was said about spherical accretion needs 
modihcation for accretion with angular momentum. 

Consider the accretion in a close binary consisting of a compact (white dwarf, 
neutron star or black hole) primary of mass Mi and a main sequence companion of 
mass M 2 . The mass ratio is dehned as q = M 2 /Mi (note: q is just as often dehned 
the other way around). 

If Ml and M 2 orbit each other in a circular orbit and their separation is a, the 
orbital frequency SI is 

Sl^ = G(Mi -bM2)/a®. 

The accretion process is most easily described in a coordinate frame that corotates 
with this orbit, and with its origin in the center of mass. Matter that is stationary 
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in this frame experiences an effective potential, the Roche potential (Ch. 4 in Frank, 
King and Raine, 1992), given by 




GM GM 1 , , 

- 

ri r2 2 


( 7 ) 


where ri 2 are the distances of point r to stars 1,2. Matter that does not corotate 
experiences very different forces (due to the coriolis force). The Roche potential 
is therefore useful only in a rather limited sense; for non corotating gas intuition 
based on the Roche geometry is usually confusing. Keeping in mind this limitation, 
consider the equipotential surfaces of (7). The surfaces of stars Mi_ 2 , assumed to 
corotate with the orbit, are equipotential surfaces of (7). Near the centers of mass 
(at low values of ())r) they are unaffected by the other star, at higher $ they are 
distorted and at a critical value <I>c the two parts of the surface touch. This is the 
critical Roche surface Sc whose two parts are called the Roche lobes. Binaries lose 
angular momentum through gravitational radiation and a magnetic wind from the 
secondary (if it has a convective envelope). Through this loss the separation between 
the components decreases and both Roche lobes decrease in size. Mass transfer starts 
when M 2 hlls its Roche lobe, and is maintained by the angular momentum loss from 
the system. A stream of gas then flows through the point of contact of the two parts 
of Sc, the inner Lagrange point Li. If the force acting on it were derivable entirely 
from (7) the gas would just fall in radially onto Mi. As soon as it moves however, 
it does not corotate any more and its orbit under the influence of the coriolis force 
is different (Fig. 1). 



Fig. 1. Roche geometry for q = 0.2, with free particle orbit from Li. Dashed: circularization 
radius. 

Since the gas at Li is very cold compared with the virial temperature, its sound 
speed is small compared with the velocity it gets after only a small distance from Li. 
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The flow into the Roche lobe of Mi is therefore highly supersonic. Such hypersonic 
flow is essentially ballistic, that is, the stream flows along the path of free particles. 

Though the gas stream on the whole follows an orbit close to that of a free 
particle, a strong shock develops at the point where the orbit intersects itself. [In 
practice shocks already develop shortly after passing the pericenter at Mi , when the 
gas is decelerated again. Supersonic flows that are decelerated by whatever means in 
general develop shocks (e.g. Courant and Friedrichs 1948, Massey, 1968). The effect 
can be seen in action in the movie published in Rozyczka and Spruit, 1993]. After 
this, the gas settles into a ring, into which the stream continues to feed mass. If 
the mass ratio q is not too small this ring forms fairly close to Mi. An approximate 
value for its radius is found by noting that near Mi the tidal force due to the 
secondary is small, so that the angular momentum of the gas with respect to Mi is 
approximately conserved. If the gas continues to conserve angular momentum while 
dissipating energy, it settles into the minimum energy orbit with the specific angular 
momentum j of the incoming stream. The radius of this orbit, the circularization 
radius r^ is determined from 

(GMirc)'/' =i. 

The value of j is found by a simple integration of the orbit starting at Li and 
measuring j at some point near pericenter. In units of the orbital separation a, r^ 
and the distance rm from Mi to Li are functions of the mass ratio only. As an 
example for q = 0.2, rm ~ 0.66a and the circularization radius Tc ~ 0.16a. In 
practice the ring forms slightly outside Tc because there is some angular momentum 
redistribution in the shocks that form at the impact of the stream on the ring. 

The evolution of the ring depends critically on nature and strength of the angular 
momentum transport processes. If sufficient ‘viscosity’ is present it spreads inward 
and outward into a disk. 

At the point of impact of the stream on the disk the energy dissipated is a 
significant fraction of the orbital kinetic energy, hence the gas heats up to a significant 
fraction of the virial temperature. For a typical system with Mi = IMq, M 2 = 
0.2Mq having an orbital period of 2 hrs, the observed size of the disk (e.g. Wood et 
al. 1989b, Rutten et ah, 1992) r^/a Ri 0.3, the orbital velocity at rj about 900 km/s, 
the virial temperature at rj is 10®K. The actual temperatures at the impact point 
are much lower, due to rapid cooling of the shocked gas. Nevertheless the impact 
gives rise to a prominent ‘hot spot’ in many systems, and an overall heating of the 
outermost part of the disk. 


4. Thin disks: properties 


4.1. Flow in a cool disk is supersonic 

The equation of motion, ignoring viscosity, in the potential of a point mass is 


(9v 

'm 


^ 1„„ GM ^ 

V • Vv = - VP - r. 


( 8 ) 


where f is a unit vector in the spherical radial direction r. To compare order of 
magnitudes of the terms, choose a position rg in the disk, and choose as typical time 
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and velocity scales the orbital time scale SIq ^ and velocity SIqCo- 

Assuming for simplicity an isothermal gas, the pressure gradient term is 

1 Tl 

-VP = —TVlnp. 

P P 

In terms of the dimensionless quantities 

f=rlro, v = v/{Qoro), 

i = Qgi, V = roV, 

the equation of motion is then 

- T - I 

All terms and quantities in this equation are of order unity by the assumptions 
made, except the pressure gradient term which has the coefhcient T /Tyir- If cooling is 
important, so that T/T^ij- ^ 1, the pressure term is negligible to Rrst approximation, 
and vice versa. Equivalent statements are also that the gas moves hypersonically on 
nearly Keplerian orbits, and that the disk is thin, as is shown next. 


4.2. Disk thickness 

The thickness of the disk is found by considering its equilibrium in the direction 
perpendicular to the disk plane. In an axisymmetric disk, using cylindrical coor¬ 
dinates (n7,(/),z), measure the forces at a point ro (n7,(/),0) in the midplane, in a 
frame rotating with the Kepler rate SIq at that point. The gravitational acceleration 
—GM/r^ f balances the centrifugal acceleration at this point, but not at some 
distance z above it because gravity and centrifugal acceleration work in different 
directions. Expanding both accelerations near ro, one finds a residual acceleration 
toward the midplane of magnitude 

9z — —^qZ. 

Assuming again an isothermal gas, the condition for equilibrium in the z direction 
under this acceleration yields a hydrostatic density distribution 

p = po{w)ex^ (“^) • 


H, the scale heiqht of the disk, is given in terms of the isothermal sound speed 
c, = by 

77 — Cs /^o • 

We define b = H/r, the aspect ratio of the disk, and find that it can be expressed in 
several equivalent ways: 


6 


r 


Sir 


= M" 



1/2 


where M is the Mach number of the orbital motion. 
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4.3. Viscous spreading 

The shear flow between neighboring Kepler orbits in the disk causes friction due 
to viscosity. The frictional torque is equivalent to exchange of angular momentum 
between these orbits. But since the orbits are close to Keplerian, a change in angular 
momentum of a ring of gas also means it must move in position. If the angular 
momentum is increased, the ring moves to a larger radius. In a thin disk angular 
momentum transport (more precisely a nonzero divergence of the angular momentum 
flux) therefore automatically implies redistribution of mass in the disk. 

A simple example (Liist 1952, see also Lynden-Bell and Pringle 1974) is a narrow 
ring of gas at some distance pq. If at t = 0 this ring is released to evolve under the 
viscous torques, one finds that it first spreads into an asymmetric hump with a long 
tail to large distances. As t ^ oo one finds that the hump flattens in such a way 
that almost all the mass of the ring is accreted onto the center, while a vanishingly 
small fraction of the gas carries almost all the angular momentum to infinity. As 
a result of this asymmetric behavior essentially all the mass of a disk can accrete, 
even if there is no external torque to remove the angular momentum. 

4.4. Observations oe disk viscosity 

Evidence for the strength of the angular momentum transport processes in disks 
comes from observations of variability time scales. This evidence is not good enough 
to determine whether the processes really have the same effect as a viscosity, but if 
this is assumed, estimates can be made of the magnitude of the viscosity. 

Cataclysmic Variables give the most detailed information. These are binaries 
with white dwarf (WD) primaries and (usually) main sequence companions (for re¬ 
views see Meyer-Hofmeister and Ritter 1993, Cordova 1993). A subclass of these 
systems, the Dwarf Novae, show semiregular outbursts. In the currently most de¬ 
veloped theory, these outbursts are due to an instability in the disk (Meyer 1990 
and references therein). The outbursts are episodes of enhanced mass transfer of the 
disk onto the primary, involving a significant part of the whole disk. The decay time 
of the burst is thus a measure of the viscous time scale of the disk (the quantitative 
details depend on the model, see Cannizzo et al. 1988): 

Cisc — ^d/^7 


where rj is the size of the disk. With decay times in the order of days, this yields vis¬ 
cosities of the order 10^® cm^/s, about 14 orders of magnitude above the microscopic 
viscosity of the gas. 

Other evidence comes from the inferred time scale on which disks around proto¬ 
stars disappear, which is of the order of 10^ years (Strom et al, 1993). 

4.5. a-PARAMETRIZATION 

The process responsible for such a large viscosity has not been identified yet. Many 
processes have been proposed, some of which demonstrably work, though often not 
with an efficiency as high as the observations of CV outbursts seem to indicate. For 
other processes, such as magnetic turbulence, no theories with sufficient predictive 
power are available. In order to compare the viscosities in disks under different 
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conditions, one introduces a dimensionless vsicosity a: 




( 10 ) 


where Cg is the isothermal sound speed as before. The quantity a was introduced 
by Shakura and Sunyaev (1973), as a way of parametrizing our ignorance of the 
angular momentum transport process (their dehnition is based on a different formula 
however, and differs by a constant of order unity). 

How large can the value of a be, on theoretical grounds? As a simple model, let’s 
assume that the shear flow between Kepler orbits is unstable to the same kind of 
shear instabilities found for flows in tubes, channels, near walls and in jets. These 
instabilities occur so ubiquitously that the fluid mechanics community considers 
them a natural and automatic consequence (e.g. DiPrima and Swinney 1981, pl44 
2nd paragraph) of a high Reynolds number: 

Re=^ 

V 


where L and V are characteristic length and velocity scales of the flow. If this 
number exceeds about 1000 (for some forms of instability much less), instability and 
turbulence are generally observed. It has been argued (e.g. Zel’dovich 1981) that 
for this reason hydrodynamic turbulence is the cause of disk viscosity. Let’s look at 
the consequences of this assumption. If due to shear instability an eddy of radial 
length scale I develops, it will rotate at a rate given by the rate of shear, < 7 , in the 
flow, here 


<7 = 


dn 



The velocity amplitude of the eddy is R = al, and a held of such eddies produces a 
turbulent viscosity of the order (leaving out numerical factors of order unity) 


I^turb — ^ 


( 11 ) 


In compressible flows, there is maximum to the size of the eddy set by causality 
considerations. The force that allows an instability to form an overturning eddy is 
the pressure, which transports information about the flow at the sound speed. The 
eddies formed by a shear instability can therefore not move faster than Cg, hence 
their size does not exceed Cg/cr Ri H. At the same time, the largest eddies formed 
also have the largest contribution to the turbulent viscosity. Thus we should expect 
that the turbulent viscosity is given by eddies with size of the order H: 

V ~ 


or 

a ~ 1. 

Does hydrodynamical turbulence along these lines exist in disks? This is currently 
not clear; some arguments are discussed briefly in section 8. 
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5. Thin Disks: equations 

Consider a thin (= cool, nearly Keplerian, cf. section 4.2) disk, axisymmetric but 
not stationary. Using cylindrical coordinates (r, (note that we have changed 

notation from w to r compared with section 4.2) we dehne the surface density T, of 
the disk as 

/ OO 

pdz^2HoPo, ( 12 ) 

- OO 

where po, Hq are the density and scaleheight at the midplane. The approximate sign 
is used to indicate that the coefficient in front of H in the last expression actually 
depends on details of the vertical structure of the disk. Conservation of mass, in 
terms of E is given by 

|(rE) + |:(rEPr) = 0. (13) 

(derived by integration the continuity equation over z). Since the disk is axisym¬ 
metric and nearly Keplerian, the radial equation of motion reduces to 


vl = GM/r. (14) 

By integrating the ())-equation of motion over height and using (12), one gets an 
equation for the angular momentum balance: 

^ ^ (15) 

where SI = v^fr, and 

/ OO 

pizdz Ri Eu. (16) 

- OO 

The second approximate equality in (16) holds if u can be considered independent 
of z. The mass flux M at any point in the disk is given by 


M = —27rrEt;r 


67rr^/^-^(uEr^/^). 

or 


(17) 


The term on the right hand side in (15) is the viscous torque. It is derived most 
easily with a physical argument (Pringle, 1981)^. 


^ If you prefer a more formal derivation, the fastest way is to consult Landau and Lifshitz (1959) 
chapter 15 (hereafter LL). Noting that the divergence of the flow vanishes for a thin axisymmetric 
disk, the viscous stress a becomes (LL eq. 15.3) 


^ik — V 


dvi dvk 
dx]^ dxi 


where 77 = pv. The viscous force is 


_ da.h dri ( dv, dvk\ 

Fr = —— = —- - - + —- + r|Av^. 

dx}^ dx}^ \dx}^ dxi ) 

The two terms on the right hand side can then be written in cylindrical or spherical coordinates 
using LL eqs. (15.15-15.18). Finally the viscous torque is computed from the ^-component of the 
force, and then integrated over 2 :. 
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Assuming that v can be taken constant with height (as long as we are not sure 
what causes the viscosity this is a reasonable simplihcation) Eqs. (12-15) can be 
combined into a single equation for E: 

This is the standard form of the thin disk diffusion equation. An important conclu¬ 
sion from this equation is: in the thin disk limit, all the physics which determines the 
time dependent behavior of the disk enters through one quantitity only, the viscosity 
10 . This is the main attraction if the thin disk approximation. 


5.1. Steady Thin Disks 


In a steady disk (d/dt = 0) the mass flux M is constant through the disk and equal 
to the accretion rate onto the central object. From (17) we get the surface density 
distribution: 


loT. = 




(19) 


where n is the inner radius of the disk and /3 is a free parameter appearing through 
the integration constant. It is related to the flux of angular momentum Fj through 
the disk: 


Fj = -Mfdn.rl 


( 20 ) 


where Sli is the Kepler rotation rate at the inner edge of the disk. If the disk 
accretes onto an object with a rotation rate SI* less than Sli, one Ends (Shakura 
and Sunyaev, 1973, Lynden-Bell and Pringle, 1974) that /3 = 1, independent of SI*. 
The angular momentum flux (torque on the accreting star) is inward (spin-up) and 
equal to the accretion rate times the specific angular momentum at the inner edge 
of the disk (eq. 20). For stars rotating near their maximum rate (SI* Ri Sli) and for 
accretion onto magnetospheres, which can rotate faster than the disk, the situation 
is different (Popham and Narayan 1991, Paczyhski 1991, Bisnovatyi-Kogan 1993, 
1994). Accreting magnetospheres, for example, can spin down by interaction with 
the disk, which in that case has a surface density distribution (19) with /3 < 1 (see 
also Spruit and Taam, 1993). 


5.2. Disk Temperature 

In this section I assume accretion onto not-too-rapidly rotating objects so that fd = 1. 
The surface temperature of the disk, which determines how much energy it loses by 
radiation, is governed primarily by the energy dissipation rate in the disk, which in 
turn is given by the accretion rate. From the first law of thermodynamics we have 

1 C' 

pT—=-divF + Ov, (21) 


where S the entropy per unit mass, F the heat flux (including radiation and any 
form of‘turbulent’ heat transport), and Qv the viscous dissipation rate. For changes 
which happen on time scales longer than the dynamical time the left hand side 
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is small compared with the terms on the right hand side. Integrating over z, the 
divergence turns into a surface term and we get 

poo 

a,T^= / Ovdz, (22) 

Jo 


where is Stefan-Boltzmann’s radiation constant = arc/4. Thus the energy 
balance is local (for such slow changes): what is generated by viscous dissipation 
inside the disk at any radius r is also radiated away from the surface at that position. 
The viscous dissipation rate is equal to = aijdvi/dxj, where aij is the viscous 
stress tensor (see footnote in section 5), and this works out to be Qy = 'J/AQp'vp. 
Eq. (22), using (19) then gives the surface temperature in terms of the accretion 
rate: 


(TyT^ = = 


GM 3M 
Stt 



(23) 


This shows that the surface temperature of a steady disk, at a given distance r from 
the accreter, depends only on the product MM, but not on the highly uncertain 
value of the viscosity. For r ^ r, we have 


% 


^- 3/4 


(24) 


These considerations only tells us about the surface temperature. The internal 
temperature in the disk is quite different, and depends on the mechanism transport¬ 
ing energy to the surface. Because it is the internal temperature that determines 
the disk thickness H (and probably also the viscosity), this transport needs to be 
considered in some detail for realistic disk models. This involves a calculation of the 
vertical structure of the disk. Because of the local (in r) nature of the balance be¬ 
tween dissipation and energy loss, such calculations can be done as a grid of models 
in r, without having to deal with exchange of energy between neighboring models. 
Schemes borrowed from stellar structure computations are used (e.g. Meyer and 
Meyer-Hofmeister 1982, Pringle et al. 1986, Cannizzo et al. 1988). 

An approximation to the temperature in the disk can be found when a number of 
additional assumption is made. As in stellar interiors radiative transport dominates 
at high temperatures. The internal temperature is then found from the Eddington 
approximation for a radiative atmosphere: 

ayT^ = -^TFy, (r > 1) (25) 

where r is the optical depth and Fy the heat flux. Assuming that most of heat is 
generated near the midplane (which is the case if v is constant with height), Fy is 
constant with height and equal to UtTs ■ Finally we assume that the opacity k is 
constant with z. The optical depth is then r = kE/ 2. With (23) this yields for the 
temperature at the midplane 

E 

64 " 


(26) 
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Assuming equation of state (1), valid when radiation pressure is small, we find with 
(19) for the disk thickness: 


n / \i/io 

= (fMieyfY (Pr C P) (27) 


where re = r/(10® cm), Mie = M/(10^®g/s), and 

/ = 1 - {n/r'Y^^ . 

From this we conclude that: i) the disk is thin, p[/r < 10“^ in X-Ray Binaries, ii) 
the disk thickness is relatively insensitive to the parameters, especially a, k and r. 
It must be stressed however that this depends fairly strongly on the assumption that 
the energy is dissipated in the disk interior. If the dissipation takes place close to 
the surface [such as in some magnetic reconnection models (Field and Rogers 1993 
and references therein)], the disk temperature will be much closer to the surface 
temperature and p[ would be even smaller in such disks. 


5.3. Radiation pressure dominated disks 

In the inner regions of disks in XRB, the radiation pressure can dominate over the 
gas pressure, which results in a different expression for the disk thickness. The total 
pressure P is 

P = P, + Pg = yT^ + Pg. ( 28 ) 

The isothermal sound speed is given by = P/p, and the relation Cg = SliJ still 
holds. For Pj- ^ Pg we get from (25), with (23) 


cH'^ 


— —KfM. 

IbTT Po 


With the rather approximate relation X = 277 po we get 


7?* 


3 K 
Sit cR^, 


fM 


3 M 
2-'Me’ 


(29) 


where 7?* is the stellar radius and Me the Eddington rate for this radius. R follows 
that the disk becomes thick near the star, if the accretion rate is near Eddington 
(though this is mitigated somewhat by the decrease of the factor /). Accretion 
near the Eddington limit is geometrically not thin any more, and in addition other 
processes such as angular momentum loss by ‘photon drag’ (Lamb, this volume) 
have to be taken into account (see also Miller, this volume). 


6. Comparison with CV observations 

The number of meaningful quantitative tests between theory of disks and observa¬ 
tions is somewhat limited since in the absence of a theory for u, it is a bit meagre 
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on predictive power. The most detailed information perhaps comes from modeling 
of CV outbursts. 

Two simple tests are possible (nearly) independently of u. These are the pre¬ 
diction that the disk is geometrically quite thin (eq. 27) and the prediction that 
the surface temperature Tg ~ r~^l^ in a steady disk. The latter can be tested in a 
subclass of the CV’s that do not show outbursts, the nova-like systems, which are 
believed to be approximately steady accreters. If the system is also eclipsing, eclipse 
mapping techniques can be used to derive the brightness distribution with r in the 
disk (Horne, 1985, 1993). If this is done in a number of colors so that bolometric 
corrections can be made, the results (e.g. Rutten et al. 1992) show in general a fair 
agreement with the r~^l^ prediction. Two deviations occur: i) a few systems show 
signihcantly flatter distributions than predicted, and ii) most systems show a ‘hump’ 
near the outer edge of the disk. The latter deviation is easily explained, since we 
have not taken into account that the impact of the stream heats the outer edge of 
the disk. Though not important for the total light from the disk, it is an important 
local contribution near the edge. 

Eclipse mappingof Dwarf Novae in quiescence give a quite different picture. Here, 
the inferred surface temperature is often nearly flat (e.g. Wood et al. 1989a). This 
is understandable however since in quiescence the mass flux depends strongly on r. 
In the inner parts of the disk it is small, near the outer edge it is close to its average 
value. With eq. (23), this yields a flatter 71 (r). The lack of light from the inner 
disk is compensated during the outburst, when the accretion rate in the inner disk is 
higher than average (see Mineshige and Wood 1989 for a more detailed comparison). 
The effect is also seen in the 2-dimensional hydrodynamic simulations of accretion 
in a binary by Rozyczka and Spruit (1993). These simulations show an outburst 
during which the accretion in the inner disk is enhanced, between two episodes in 
which mass accumulates in the outer disk. The bolometric surface brightness in this 
simulation is shown in Fig. 2. 


7. Comparison with LMXB observations 

In low mass X-ray binaries a complication arises because of the much higher lumi¬ 
nosity of the accreting object. Since a neutron star is roughly 100 times smaller than 
a white dwarf, it produces 100 times more luminosity for a given accretion rate. 

Irradiation of the disk by the central source leads to a different surface temper¬ 
ature than predicted by (23). The central source (star plus inner disk) radiates the 
total accretion luminosity GMM/i?* (assuming sub-Eddington accretion, see sec¬ 
tion 2). If the disk is concave, it will intercept some of this luminosity. If the central 
source is approximated as a point source, and the angle between the disk surface 
and the radial direction is e, the irradiating flux on the disk surface is 

c = diJ/dr — H/r, (30) 
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Fig. 2. Bolometric surface brightness in numerical simulations of Rozyckzka and Spruit, 1993. 
Long dashed: during, solid: before, short dashed: after the outburst. Note hump at the disk edge, 
due to the impact of the stream. 


where is the flux generated internally in the disk, given by (23). On average, the 
angle e is of the order of the aspect ratio 6. With / Ri 1, and our fiducial value 6 Ri 
5 10“^, we And that irradiation dominates for r > 10®cm. This is compatible with 
observations (for reviews see Mason 1990, Van Paradijs and McClintock 1993), which 
show that the optical and UV are completely dominated by reprocessed radiation. 
When irradiation by an external source is included in the thin disk model, the energy 
balance equation (23) becomes 




GM?>M 

Stt 


/ + (1 


^).^irr7 


(31) 


where a is the X-ray albedo of the surface, i.e. 1 — a is the fraction of the incident flux 
that is absorbed in the optically thick layers of the disk (photons absorbed higher 
up only serve to heat up the corona of the disk). In the region where the irradiating 
flux dominates the first term on the rhs can be neglected in computing the surface 
temperature. For Tirr we need the disk thickness FI. It turns out that the irradiating 
flux, even when it is larger than the internally generated heat flux, has only a small 
effect on the disk thickness, as long as Firr/Fd < t, where r is the optical thickness 
of the disk (Lyutyi and Sunyaev 1976, Vrtilek et al. 1991). Here, we are assuming 
again that the energy is transported vertically by radiation. The reason for this is 
the same as in radiative envelopes of stars, which are also insensitive to the surface 
boundary condition. In the reprocessing region of the disks of LMXB, the conditions 
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are such that _fd << -Flrr < tF^, hence we must use eq. (27) for H. This yields 
e = (21/20)-ff/r Ri 5 10“^, and Tg ~ ®, and we still expect the disk to remain thin. 

From the paucity of sources in which the central source is eclipsed by the com¬ 
panion one deduces that the companion is barely or not at all visible from the inner 
disk, presumably because the outer parts of the disk are much thicker than expected 
from the above arguments. This is consistent with the observation that the charac¬ 
teristic modulation of the optical light curve due to irradiation of the secondary’s 
surface by the X-rays is not very strong in LMXB (with the exception of Her X-1, 
which has a large companion). The place of the eclipsing systems is taken by the 
so-called ‘Accretion Disk Corona’ (ADC) systems, where shallow eclipses of a usu¬ 
ally extended X-ray source are seen (for reviews of the observations, see Lewin et al. 
1995). The conclusion is that there is an extended X-ray scattering ‘corona’ above 
the disk. It scatters a few per cent of the X-ray luminosity. 

What causes this corona and the large inferred thickness of the disk is not clear at 
present. The thickness expected from disk theory is a rather stable number. To ‘sus¬ 
pend’ matter at the inferred height of the disk forces are needed that are much larger 
than the pressure forces available in an optically thick disk. A thermally driven wind, 
produced by X-ray heating of the disk surface, has been invoked (Begelman et ah, 
1983, see also Horn and Meyer, 1993). For other explanations, see van Paradijs and 
McClintock (1995). Perhaps a magnetically driven wind from the disk, such as seen 
in protostellar objects (e.g. Kbnigl and Ruden 1993) can explain both the shielding 
of the companion and the scattering. Such a model would resemble magnetically 
driven wind models for the broad-line region in AGN (e.g. Emmering et ah, 1992). 

7.1. Transients 

Soft X-ray transients (also called X-ray Novae) are believed to be binaries similar to 
the other LMXB, but somehow the accretion is episodic, with very large outbursts 
recurring on time scales of decades (sometimes years). There are many black hole 
candidates among these transients (see van Paradijs and McClintock 1995 for a 
review). As with the Dwarf Novae, the time dependence of the accretion in transients 
can in principle be exploited to derive information on the disk viscosity, assuming 
that the outburst is caused by an instability in the disk. The closest relatives 
of soft transients among the White Dwarf plus main sequence star systems are 
probably the WZ Sge stars (van Paradijs and Verbunt 1984, Bailyn 1992), which 
show (in the optical) similar outbursts with similar recurrence times (cf. Warner 
1987, O’Donoghue et ah 1991). Like the soft transients, they have low mass ratios 
(q < 0.1). For a given angular momentum loss, systems with low mass ratios have 
low mass transfer rates, so the speculation is that the peculiar behavior of these 
systems is somehow connected with a low mean accretion rate. 

7.2. Disk Instability 

The most developed model for outbursts is the disk instability model of Osaki (1974), 
Hoshi (1979), Smak (1984), Meyer and Meyer-Hofmeister (1981), Faulkner et ah 
(1983) (for recent discussions see King 1995, Osaki 1993). In this model the insta- 
blity that gives rise to cyclic accretion is due to a temperature dependence of the 
viscous stress. In any local process that causes an effective viscosity, the resulting 
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a- parameter will be a function of the main dimensionless parameter of the disk, 
the aspect ratio H/r. If this is a sufficiently rapidly increasing function, such that 
a is large in hot disks and low in cool disks, an instability results by the following 
mechanism. Suppose we start the disk in a stationary state at the mean accretion 
rate. If this state is perturbed by a small temperature increase, a goes up, and by 
the increased viscous stress the mass flux M increases. By (23) this increases the 
disk temperature further, resulting in a runaway to a hot state. Since M is larger 
than the average, the disk empties partly, reducing the surface density and the cen¬ 
tral temperature (eq. 26). A cooling front then transforms the disk to a cool state 
with an accretion rate below the mean. The disk in this model switches back and 
forth between hot and cool states. By adjusting a in the hot and cool states, or by 
adjusting the functional dependence of a on H/r, outbursts are obtained that agree 
reasonably with the observations of soft transients (Lin and Taam 1984, Mineshige 
and Wheeler, 1989). A rather strong increase of a with H/r is needed to get the 
observed long recurrence times. 

Another possible mechanism for instability has been found in 2-D numerical 
simulations of accretion disks (Blaes and Hawley 1988, Rozyczka and Spruit 1993). 
The outer edge of a disk is found, in these simulations, to become dynamically 
unstable to a oscillation which grows into a strong excentric perturbation (a crescent 
shaped density enhancement which rotates at the local orbital period). Shock waves 
generated by this perturbation spread mass over most of the Roche lobe, at the same 
time the accretion rate onto the central object is strongly enhanced. This process is 
different from the Osaki and Hoshi mechanism, since it requires 2 dimensions, and 
does not depend on the viscosity (instead, the internal dynamics in this instability 
generates the effective viscosity that causes a burst of accretion). 

7.3. Other Instabilities 

Instability to heating/cooling of the disk can be the due to several effects. The 
cooling rate of the disk, if it depends on temperature in an appropriate way, can 
cause a thermal instability like that in the interstellar medium. Other instabilities 
may result from the dependence of viscosity on conditions in the disk. For a general 
treatment see Piran (1978), for a shorter discussion see Treves et ah, 1988. 


8. Sources of Viscosity 

The high Reynolds number of the flow in accretion disks (of the order 10^^ in the 
outer parts of a CV disk) would, to most fluid dynamicists, seem an amply sufficient 
condition for the occurrence of hydrodynamic turbulence. A theoretical argument 
against such turbulence often used in astrophysics (Kippenhahn and Thomas 1981, 
Pringle 1981) is that in cool disks the gas moves almost on Kepler orbits, which are 
quite stable (except the orbits that get close to the companion). This stability is re¬ 
lated to the known stabilizing effect that rotation has on hydrodynamical turbulence 
(Bradshaw 1969, for a recent discussion see Tritton 1992). Kippenhahn and Thomas 
also point out that the one laboratory experiment that comes close to the situation 
in accretion disks, namely the rotating Couette flow, does not become unstable for 
parameters like in disks (for the rather limited range in Reynolds numbers available). 
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A (not very strong) observational argument is that hydrodynamical turbulence as 
described above would produce an a that does not depend on the nature of the disk, 
so that all objects should have the same value. This is unlikely to be the case. From 
the modeling of CV outbursts one knows, for example, that a probably increases 
with temperature (more accurately, with H/r, see previous section). Also, there are 
indications from the inferred life times and sizes of protostellar disks (Strom et al. 
1993) that a may be rather small there, ~ 10“^, whereas in outbursts of CV’s one 
infers values of the order 0.1 — 1. 

The indeterminate status of the hydrodynamic turbulence issue is one of the 
most annoying problems in disk theory. Numerical simulations may be of some 
help, even though Reynolds numbers as high as those in disks can not be reached. 
The astrophysical approach has been to circumvent the problem by Rnding plausible 
alternative mechanisms that might work just as well. 

Among the processes that have been proposed repeatedly as sources of viscosity 
is convection due to a vertical entropy gradient (for recent work see Kley et al. 
1993), which may have some limited effect in convective parts of disks. Another 
class are waves of various kinds. Their effect can be global, that is, not reducible 
to a local viscous term because by traveling across the disk they can communicate 
torques over large distances. For example, waves set up at the outer edge of the disk 
by tidal forces can travel inward and by dissipating there can effectively transport 
angular momentum outward (e.g. Narayan et al. 1987, Spruit et al. 1987). A 
nonlinear version of this idea are selfsimilar spiral shocks, observed in numerical 
simulations (Sawadaet al. 1987) and studied analytically (Spruit 1987). Such shocks 
can produce accretion at an effective a of 0.01 in hot disks, but are probably not very 
effective in disks as cool as those in CV’s and XRB. A second non-local mechanism 
is provided by a magnetically accelerated wind originating from the disk surface 
(Blandford 1976, Bisnovatyi-Kogan and Ruzmaikin 1976, Lovelace 1976, Blandford 
and Payne 1982, for recent reviews see the papers in Lynden-Bell 1993). In principle, 
such winds can take care of all the angular momentum loss needed to make accretion 
possible in the absence of a viscosity (Kbnigl, 1989). The attraction of this idea is 
that magnetic winds are a strong contender for explaining the strong outflows and 
jets seen in protostellar objects and AGN. It is not yet clear however if, even in 
these objects, the wind is actually the main source of angular momentum loss. In 
sufhciently cool or massive disks, selfgravitating instabilies of the disk matter can 
produce internal friction. Paczynski (1978) has proposed that the resulting heating 
would limit the instability and keep the disk in a well defined moderately unstable 
state. The angular momentum transport in such a disk has been modeled by several 
authors (e.g. Lin and Pringle, 1990). Disks in XRB are always too hot for selfgravity 
to play a role. 

Magnetic forces can be very effective at transporting angular momentum. If it 
can be shown that the shear flow in the disk produces some kind of small scale 
fast dynamo process, that is, some form of magnetic turbulence, a large effective 
a may be expected. An elegant demonstration that disks do indeed generate such 
fields on a dynamical time scale has recently (Balbus and Hawley 1991, 1992) come 
with the revival of an old instability meachanism in magnetic shear flows (Velikhov 
1959, Chandrasekhar 1961). This instability rapidly stretches a weak initial held 
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into small loops in the r — z plane, which are then stretched into an azimuthal 
held by the shear. The importance of the mechanism is that it will amplify a 
weak held even if it has initially no component that can be stretched by the shear 
how. The nonlinear development of this instability is almost certainly highly 3- 
dimensional, involving a mostly small-scale magnetic held. It is not clear that the 
instability still plays a signihcant role in this highly dynamic situation; the net 
effect could well be the magnetic turbulence whose properties have been guessed 
previously by many authors (e.g. Eardley and Lightman 1975, Pudritz 1981, Meyer 
and Meyer-Hofmeister 1982). The main difficulty with these models has been the 
lack of reliability of the end result, the effective a. It is conceivable that by 3-D 
numerical magnetohydrodynamic simulations one may soon get sufficient insight to 
derive useful numbers. 

9. Beyond thin disks 

Ultimately, much of the progress in developing useful models of accretion disks will 
depend on detailed numerical simulations in 2 or 3 dimensions. In the disks one 
is interested in, there is usually a large range in length scales (in LMXB disks, 
from less than the 10km neutron star radius to the more than 10®km orbital scale). 
Correspondingly, there is a large range in time scales that have to be followed. 
This not technically possible at present and in the foreseeable future. In numerical 
simulations one is therefore limited to studying in an approximate way aspects that 
are either local or of limited dynamic range in r,t (for examples, see Hawley 1991, 
Rozyczka and Spruit 1993). For this reason, there is still a need for approaches that 
relax the strict thin disk framework somewhat without resorting to full simulations. 

9.1. ‘Slim’ disks 

Due to the thin disk assumptions, the pressure gradient does not contribute to the 
support in the radial direction and the transport of heat in the radial direction is 
negligible. Some of the physics of thick disks can be included in a fairly consistent 
way in the so-called ‘slim disk’ approximation (Abramowicz et ah, 1988). It has the 
advantage of remaining one-dimensional, but allows the inclusion of the radial pres¬ 
sure gradient (so that for example sound waves can occur) and radial heat transport. 
It still ignores all multidimensional processes, but allows one to study those aspects 
of the physics of thick disks that are approximately one-dimensional. 

9.2. Boundary layers 

In order to accrete onto the star, the disk matter must dissipate an amount of energy 
given by 

^^^(l-^2*/^2k(i^*))^ (32) 

Z-tC^ 

The factor in brackets measures the kinetic energy of the matter at the inner edge 
of the disk (r = i?*), in the frame of the stellar surface. Due to this heating the disk 
inflates into a ‘belt’ at the equator of the star, of thickness H and radial extent of 
the same order. Equating the radiation emitted from the surface of this belt to (32) 
the surface temperature of the belt, assuming optically thick conditions and a 
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slowly rotating star (Sl*/Slk ^ 1) is given by 


To find the temperature inside the belt and its thickness, use eq. (25). The value 
of the surface temperature is higher, by a factor of the order than the 

simplest thin disk estimate (23, ignoring the (r/n)^/^ factor). In practice, this 
works out to a factor of a few. The surface of the belt is therefore not very hot. The 
situation is quite different if the boundary layer is not optically thick (Pringle and 
Savonije 1979). It then heats up to much higher temperatures. Analytical methods 
to obtain the boundary layer structure have been used by Regev and Hougerat 
(1988), numerical solutions of the slim disk type by Narayan and Popham (1993), 2- 
D numerical simulations by Kley (1991). These considerations are primarily relevant 
for CV disks; in accreting neutron stars, the dominant effects of radiation pressure 
have to be included. 
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